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ABSTRACT 

We measure the integrated contributions of dusty AGB stars and other luminous red sources to 
the mid-IR luminosities of 6 galaxies (M81, NGC 2403, NGC 300, M33 and the Magellanic Clouds). 
We find the dusty AGB stars contribute from 0.6% (M81) to 5.6% (SMC) of the 3.6 fim flux and 
1.0% (M81) to 10.1% (SMC) of the 4.5 /zm flux. We find a trend of decreasing AGB contribution 
with increasing galaxy metallicity, luminosity and mass and decreasing SSRF. However, these galaxy 
properties are strongly correlated in our sample and the most likely explanation of the trend is galaxy 
metallicity. The bright, red sources represent a smaller fraction of the luminosity, ~1.2% at 3.6 /im, 
however their dust is likely cooler and their contributions are likely larger at longer wavelengths. 
Excluding the SMC, the contribution from these sources correlates with the specific star formation rate 
as we would expect for massive stars. In total, after correcting for dust emission at other wavelengths, 
the dust around AGB stars radiates 0.1-0.8% of the bolometric luminosities of the galaxies. Thus, 
hot dust emission from AGB and other luminous dusty stars represent a small fraction of the total 
luminosities of the galaxies but a significant fraction of their mid-IR emissions. 

Subject headings: stars: AGB - galaxies: stellar content - galaxies: individual(M81, NGC 300, 
NGC 2403, M33) - Magellanic Clouds 



1. INTRODUCTION 

Understanding the emissions from AGB stars is im- 
portant for stellar evolution, models of galaxy spectral 
energy distributions (SEDs) and inferences about stellar 
populations and their evolution. AGB stars are H and 
He shell burning stars that evolve from 0.8-8 M Main 
Sequenc e stars. During the AGB phase, which lasts 1- 
13 Myr (Vassiliadis fc Wood|1993 ), the star can undergo 
significant mass-loss, particularly in the Thermal-Pulsing 
AGB phase (TP-AGB, M < 5M ). During this last 0.2- 
2 Myr of the AGB phase, the AGB star experiences ther- 
mal pulses from shell-shell interactions that drive high 
mass loss rate winds conducive to the formation of dust. 
The short lifetime of this phase makes these stars rel- 
atively rare, so they are generally absent from the star 
clusters used to calibrate stellar evolution models. As a 
result, this is one of the most uncertain phases of stellar 
evolution. Their large range of mass loss rates combined 
with dust formation mean that AGB stars can be impor- 
tant over a broad range of wavelengths from the optical 
to the mid-IR. 

For distant galaxies, the properties of stellar popula- 
tions must be inferred from spectra or the SED of en- 
sembles of stars. Since AGB stars represent a non-trivial 
fraction of the luminosity of ^Gyr old stellar popula- 
tions, any uncertainties in the AGB phase propagate 
into uncertainties in overall galaxy properties because 
the stellar population synthesis (SPS) model must as- 
sume some t reatment of the (TP) AGB stars. Maraston 
et al.| ( |2006[ ) shows how different calibrations of the TP- 
AGB stars m the SPS models can cause large, systematic 
differences in the fits to SEDs and the resulting estim ates 
of a galaxy's mass and age. Conroy et al. (2010) and 
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Conroy & Gunn (2010) show in a more general way how 
TP-AGB stars can effect the determination of a galaxy's 
properties. 

Much progress has been made recently in moving be- 
yond the Milky Way and Magellanic Cloud star clus- 
ters to calibrate the AGB phase of stellar evolution (e.g . 
Kriek et al.|2010[ |Meidt et al.|2012bl|Zibetti et al.|2012[ ). 



In particular, [Melbourne et al.| ( |2012 ) used N1R obser 



vations to build on the optica l study of 12 metal poor 



galaxies by Girardi et al. ( |2010 ) to constrain the lifetimes 
of AGB stars based on trie number and luminosities of 
the AGB stars. While they found that the numbers of 
observed AGB stars agreed with the updated ( |Girardi et 
al][2l)Iol) SPS models given their uncertainties, the pre 



dieted luminosities are too large. They go on to demon- 
strate how more accurate lifetimes and luminosities are 
needed t o be able to accura tely describe galaxies at high 
redshift. Boyer et al. ( 2009 ) examined AGB stars in eight 
Local Group dwarf irregular galaxies in the Spitzer IRAC 
bands. They examined the number of AGB stars, their 
mass-loss rates and how this mass returning to the ISM 
could effect the current SFR. They also suggest that op- 
tical AGB searches will miss 30%-40% of AGB stars due 
to dust obscuration. 

In this paper we focus on the contribution of hot dust 
around luminous stars, particularly AGB stars, to the 
mid-IR emission of 6 nearby galaxies of varying prop- 
erties. In addition to AGB stars, these sources include 
dusty young star clusters and other luminous stars that 
have experienced dense, episodic winds or eruptive out- 
bursts. We use archival Spitzer IRAC 3.6 and 4.5 /j,m ob- 
servations of M33, M81, NGC 300, NGC 2403, the SMC 
and the LMC to survey these populations. Section 2 de- 
scribes the data and its analysis. We are interested in the 
mid-IR emissions of dusty stars that have red [3.6] — [4.5] 
colors indicating the presence of hot circumstellar dust. 
Section 3 presents the data analysis, including the source 
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Fig. 1. — Spitzer 3.6/xm image of M33. The solid curve shows 
the region used for the galaxy and the dashed lines outline the 
background regions. 

classification and background corrections. In Section 4 
we determine the fraction of the [3.6] and [4.5] flux con- 
tributed by the AGB stars and examine the trends in 
the AGB flux fraction with galaxy properties. Our final 
summary is in Section 5. 

2. DATA AND ANALYSIS 



We used archival Spitzer IRAC ( |Fazio et al.||2004[ ) 3.6 
and 4.5 fim images of M81, NGC 300, NUC 2403 and 
M33. For the SMC and LMC we used the SAG E cat- 
alogs ( |Gordon et ~aLl|2011| |Meixner et~al] |2006|) . The 
M33 data w ere the six co-added epochs f rom McQuinn 
et al. l (|2~007 ) that were also reanalyzed by | Thompson et 
a ppOOlIlr The Local Volume Legacy (LVL) Survey ([ Dale 



etaT 2009) data were used for NGC 300. For M8Tand 
NUG1303: we used the SINGS Legacy Survey flKenni- 
cutt et aL||2003 ) data. The M33 data have a pixel scale 
of 1'.' 2 pixel" \ while NGC 300, M81 and NGC 2403 have 
a pixel scale of (X'75 pixel -1 . We initially examined the 
[5.8] and [8.0] band data as well, but resolution-induced 
confusion made this problematic and we decided to ex- 
clude these bands from further analysis. 

For each mosaic wc defined the galaxy as an ellipse 
centered at the same coordinates and with the same po- 
sition angle and ellipticity as used by the LVL Survey. 
The size of the ellipse was generally limited by the field 
of view of the observation. For M33 we used a semi- 
major axis of 13.'5. For NGC 300, M81 and NGC 2403 
we used a semimajor axis of 10'2. For the SMC and 
LMC we used circles visually centered on the galaxies 
that were 2.0° and 3.5° in radius, respectively. We draw 
the sources and determine global fluxes for the galaxies 
from these areas. Figure [T] illustrates this for M33. 

We identified sources on the IRAC mosaics using the 
DAOPHOT/A LLSTAR PSF fitting photometry package 
( Stetson[l9 92 ) . A minimum of 5 relatively isolated stars 
were used to determine the PSF, and sources were re- 
quired to be at least 3a above the background at both 
[3.6] and [4.5] with a positional match of less than 1 pixel. 
Stars with an apparent magnitude of [3.6] < 8.5 mag 
were masked as foreground stars. We chose this magni- 



tude because it would include any possible r\ Carinae in 
our nearest galaxy with images, M33, while removing the 
brightest foreground stars. We masked a 15 to 25 pixel 
radius, depending on the brightness of the masked star. 
Using the distances in Table [TJ the fluxes were trans- 
formed to an absolute scale of vL v in solar luminosities. 

Table [T] summarizes the properties of the galaxies. 
Most of these are taken from the literature, as reported 
in the table caption. The star forma tion rates (SFR) 
were derived using Ha measurements dK ennicutt et ah 



2008|) and calibrated to a SFR following |Meurer et al 
. The stellar masses for the four more distant 



galaxies were determined from the 2MASS ( pkrutskie et 
al. |2006 ) K band magnitudes measured by I Dale et al . 
2009 ) , assuming a mass-to- light r atio of 0.95 flBell et a 



Their metallicities are from Zaritsky et al. ( |1994[ ) 



and were calculated at half the semimajor axis of the 
galaxy. While there is some controversy over the actual 
galaxy metallicities and gradients, it is the differential 
metallicities that are important. For example, Urbancja 



et al. ( 2005 1 present a central metallicity and gr adient for 
NGC 3U0 th at are much lower than is found in |Zaritsky| 
et al.| (fl"994|), however when they fit their data with the 
Zaritsky et al. ( 1994 ) i? 2 3 calibration, the metallicity esti- 
mates differ by only 0.05 at half the semimajor axis of the 
galaxy. Such differences will not effect our conclusions. 
The masses for the S MC and LMC wer e determined using 
the 3.6/im flux from Dale et al. (2009), assuming a mass 
to light ratio found by taking the uncertainty weighted 
averag e of the 24 Sa galaxies from|Falc6n-Ba rroso et ah] 
fl2011|) . The metallicities of the SMC and LMC are taken 



from IPei mbert &; Torres-Peimbert (1976) and Pagel et 
al. ( 1978| , respectively. The specific star formation rate 



(5SFR) is the SFR per unit stellar mass. 

The total mid-IR luminosity of each galaxy was de- 
termined in two steps, with the measurements show in 
Table [2] We first used the IRAF Ellipse task to mea- 
sure the flux inside the previously defined elliptical re- 
gions. To correct this flux for foreground and background 
emission, we next estimated the mean surface brightness 
in a background region outside the ellipse defining the 
galaxy and within the mosaic's area of equal exposure 
time. The background areas for M33, M81, NGC 300 and 
NGC 2403 are small sections located as far as possible 
from the galaxy given the mosaic geometry (see Figure 
[l]). The backgrounds for the LMC and SMC were deter- 
mined from annuli with sections cut out where the expo- 
sure times differed from the central regions. For the LMC 
we use an annulus from 3.7° to 4.0° and for the SMC the 
annulus is from 2.3° to 2.5°. We then subtracted this 
mean background surface brightness multiplied by the 
area of the signal region from the flux found by the El- 
lipse task. We have uniformly masked the rarer bright 
foreground stars, about 10 stars per image. We cannot 
mask too many stars before starting to lose the bright- 
est mid-IR sources in the target galaxies. However, we 
tested masking a larger number of stars, and found only 
a modest effect on the luminosity estimates. Our back- 
grou nd corrected galax y luminosities are in good agree- 



ment ipaleetaL (2009). The four more distant galaxies 
agree within I c, except for the 4.5 /jm luminosity of M33 
and NGC 300, which agree within 2a. The LMC is the 
most discrepant, at ~ 4c. The difference is likely caused 
by our smaller galaxy aperture. 
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Fig. 2. — CMDs for all 6 galaxies. The AGB stars are in green, 
the luminous dusty stars are in red, and the normal stars are in 
blue. 



3. CMDS AND SOURCE CLASSIFICATION 

We classified the sources using the M4.5 and [3.6] — [4.5] 
color magnitude diagram (CMD). Figure [2] shows the 
CMDs with classification criteria for all 6 galaxies. We 
used a conservative lower luminosity limit of M4.5 < 
—8.75 mag. This limit is to ensure accurate photome- 
try and matches between bands for all galaxies. Color 
errors dominate any attempt to classify lower luminosity 
sources in the more distant galaxies. We divided the stars 
into three regions using the boundaries defined in Table 
[3l The boundaries were largely determined from the M33 
(JMD, which has a clear and well-populated dusty AGB 
sequence. We first defined a region consisting of normal 
stars with little dust (blue). In addition to the lower 
luminosity limit, we also used an upper limit to elimi- 
nate bright, foreground stars, a blue limit to eliminate 
peculiar outliers and a red limit to remove stars with hot 
dust. Next we defined an AGB region (green) based on 
the CMD of M33. Finally, we defined a region to include 
any other luminous dusty stars (red) , which we will refer 
to as the red star region. The upper luminosity limit was 
allowed to be higher here in order to include objects such 
as 77 Carinae and dusty star clusters. 

The total flux associated with each class of sources 
is then found by summing the fluxes of the individual 
stars. There will also be contamination from both fore- 
ground stars and background extragalactic sources. As 
with the total fluxes, we estimated this contamination 
using the background regions defined for each galaxy. 
We built CMDs for each background region, computed 
the mean fluxes, and subtracted them from the total flux 
after adjusting for the differences in sampling area. Ta- 
ble [4] shows the percentage of source light that is con- 
sidered background for each of our source classes in each 
galaxy. The background corrections can be quite large. 
Figure [3] shows the background region CMD for M33 that 
has been corrected for area by randomly drawing sources 
from the background region and adding a small random 
scatter to the magnitudes so that the effective area is the 
same as for the galaxy CMD in Figure [2| We can still 



Fig. 3. — The CMD of the M33 background regions. To correct 
the density of the points for the smaller area of the background 
regions we have corrected the total number of objects for the area. 
We randomly drew an area-corrected number of sources from the 
background region and applied a small random scatter to the mag- 
nitudes so that the background region CMD can be visually com- 
pared to Figure [2] 



see a weak AGB feature in this background CMD, so the 
background region still includes stars in M33. Ideally we 
would use background regions better separated from the 
galaxies to eliminate this problem, but we are limited by 
the observations and have chosen regions as far from the 
galaxies as the data allow. While this implies some over- 
subtraction, we apply this background correction to both 
the summed luminosity of the sources and the estimates 
of the total galaxy luminosity. Since we ultimately fo- 
cus on the ratio of these two quantities, the final results 
should be little effected by any oversubtraction. For ex- 
ample, the AGB luminosity fraction determined for M33 
using an elliptical background region extending from 17'.6 
to 20'.0 rather than our standard region shown in Fig- 
ureflj changes the AGB luminosity fraction by less than 
0.1% in both IRAC bands. We carried out similar tests 
for each galaxy. M33 showed the largest changes, and the 
differences are always smaller than our estimated statis- 
tical uncertainties. 

We estimate the statistical uncertainties using boot- 
strap resampling. For the total galaxy flux for M33, M81, 
NGC 300 and NGC 2403, we calculated the uncertain- 
ties in the background fluxes by bootstrap resampling the 
image in 5x5 pixel squares. Using patches comparable 
in size to point sources will roughly include the Poisson 
fluctuations of point sources in the uncertainties. For the 
LMC and SMC, where we were working from catalogs, 
we bootstrap resampled the sources both from the back- 
ground and galaxy regions. We also bootstrap resampled 
the AGB, red, and blue sources in both the galaxy and 
the background region to determine the uncertainties in 
both the total stellar flux for each category and their 
background contamination correction. 

Finally, we calculated the fraction of the total lumi- 
nosity that is attributed to each type of star for both the 
original and the background-corrected values. Table [5] 
shows the fraction of the luminosity that the AGB, rea 
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Fig. 4. — The fraction of the 3.6/im (filled triangles) and 
4.5/^m (open circles) luminosity produced by stars in the AGB 
region as a function of galaxy luminosity (top left), stellar mass 
(top right), specific star formation rate (lower left) and metallicity 
(lower right). 

and normal stars contribute to eachgalaxy without any 
background correction, while Table [6] shows these quan- 
tities corrected for background contamination. We will 
discuss only the background corrected values. 

4. RESULTS 

We start by discussing the AGB region, which accounts 
for 0.6% (M81) to 5.6% (SMC) of the 3.6 flux and 
1.0% (M81) to 10.1% (SMC) of the 4.5 /an flux. As 
shown in Figure |4j we see a general trend of decreas- 
ing AGB contribution with increasing galaxy luminosity, 
mass and metallicity and an increasing AGB light frac- 
tion with increasing SSFR. The most likely interpreta- 
tion is that the AGB fraction of the mid-IR luminosity 
is higher for the lower metallicity galaxies. While it ap- 
pears the mid-IR flux contribution also correlates with 
luminosity, mass, and SSFR, we believe this is just a con- 
sequence of the lower metallicity galaxies having higher 
SSFRs and lower masses. Since the AGB stars are 0.5-2 
Gyr old, the AGB luminosity should have no intrinsic 
correlation with the current SFR or SSFR. The star for- 
mation within that 1.5 Gyr window will tightly correlated 
with the AGB luminos ity contribution (Kelson & Holden 
2010 Kriek et al. 2010 ), but we lack complete star forma- 
tion histories tor all the galaxies in our full sample. Sim- 
ilarly, the stellar mass depends on the integrated SFR, 
which should also show no direct correlation. 

The dependence of AGB dust production and the re- 
sulting IR emission on initial metallicity is an area of 
active study, with somewhat conflicting results. It is a 
difficult problem, combining stellar evolution, chem istry, 
wind formation a nd dust properties. For example, Ven- 
tura et al. ( 2012[ ) examined models of dust production 
by AGB stars at LMC and SMC metallicities. They find 
that for the more massive stars that undergo Hot Bottom 
Burning (8M Q > M > 3.5Mq), a higher metallicity re- 
sults in more dust production. For lower mass stars that 
undergo third dredge up (3.5M© > M > 1.5M ), they 
find that the dust production is nearly independent of 



metallicity. The more massive stars become oxygen rich 
at the surface, in which case dust formation is limited 
by the availability of silicon, which is controlled by the 
metallicity. For the low mass stars, third dredge up leads 
to an enhanced surface abundance of carbon and thus 
dust formation that d epends little on initial metallicity. 
Ventura et al. (2012) in the end argue that theoretical 
predictions of dust production around lower mass AGB 
stars are not robust due to uncertainties in the amount 
of mixing and the exten t of the third dredge up. Simi 



lar results are found by|Marigo et al. (|2008 ) , Bowen fe; 



Willson (1991) and IMatsuura et al.l (2007). Groenewe 



gen et al. 1 2009| ) compare their sample of LMC and SMC 
AGB stars to models for mass loss, taking different dust 
grain composition into account, and compare to Galactic 
estimates to find that there is no strong dependence of 
mass-loss rate on metallic ity to factors o f order 2-4. This 
is broadly consistent with Wachter et al. ( 2008[ ), who find 
the mass loss rate of the Milky Way AGB stars to be a 
factor of 2 higher than that of the SMC. While higher 
metallicity may result in more dust formation, this dust 
is dominated by silicates which have a lower opacity than 
graphitic dust. These issues make it difficult to translate 
these trends into obser vational predictions. 

The Maraston ( 2005 1 models of the AGB c ontribution 
are based on the fuel consumption theorem (Renzini & 



Voli 1981 ) and assume mass loss is not significantly el 
tected by metallicity. Instead, it is the surface abun- 
dances that are effected, so that the number ratio of 
oxygen-rich to carbon-rich AGB stars depends on metal- 
licity. A metal-poor population is expected to have more 
carbon-rich stars. A metal-poor AGB star has a lower 
abundance of oxygen in its envelope, so less carbon has 
to be dredged up before the oxygen has all been bound 
in CO and the residual carbon can gp on to form dust 
( Iben fc Renzini|198 3 Mara ston et al.|2006 ). In this sce- 
nario, a metal-poor population could resu l t in a higher 



mid-IR luminosity. Bird & Pinsonneault (2011) exam 



ined the energy output of AGB stars using fuel consump 
tion theory, but find insufficient calibration data to de 
termine an accu rate metallicity dependence. The earlier 
Pado va models (|Marigo & GirardiJ[2007 



Marigo et al. 



2008 1 have more carbon-rich stars at low metallicity due 



to the higher efficiency of third dredge up, while the more 



recent |Girardi et al. (2010) models have increased mass 
loss rates and shorter AGB phase lifetimes at low metal- 
licity. While the effect of metallicity is not fully explored 
by these studies, the increased dust production and mass 
loss rate suggested at low metallicity would result in a 
higher mid-IR contribution, as we find in our sample. 

Previous observational constraints on the metallicity 
depende nce of the dust produ ction also have mixed con- 



clusions. Meidt et al. (2012b) found that lower metallic 



ity clusters in M 100 have more dusty AGB stars, arguing 
that since the dust optical depth depends upon whether 
an AGB star is O-rich or C-rich, the AGB contribution 
to the cluster depends upon metallicity. They also find 
larger mass-loss rates for younger clusters. Since the 
younger clusters also have a higher metallicity, they can- 
not address the overall effect of metallicity on the AGB 
contribution. Looking at the AGB population of eight 



local group dwarf irregular galaxies, Boyer et al. (2009) 
find that all galaxies in their sample have obscured AGB 
stars. For AGB stars with optical counterparts, they 
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find that the more metal-rich galaxies have a larger pop- 
ulation of red, dust enshrouded AGB stars. However, 
requiring optical counterparts biases the sample against 
dustier stars. An alternate explanation for the increased 
number of dusty stars in the more metal-rich galaxies is 
the age of the populations, because the higher metallic- 
ity galaxies in their sample also have more recent (< 3 
Gyr) star formation. They find no clear correlation be- 
tween metallicity and optical completeness, as would be 
expected if metallicity were the dominating factor in the 
d egree of obscuration. 

|Meidt et al.| ( |2012a[ ) examined the contribution from 
red and AGB stars to the luminosity of M81. They 
find dust contributes 5.0±0.9% and a smaller contribu- 
tion from red and AGB stars, at only 0.10±0.02% of the 
3.6 /jm luminosity, rather than our significan tly higher 
estimate of 0.58±0.02% for the AGB stars. |Meidt et| 



al. (2012a) used a very different method to estimate the 
percentage of the light due to AGB and intermediate 
age stars. Essentially, they divide the 3.6 and 4.5 /im 
emission into that from stars and that due to "contami- 
nating" sources that contribute to the mid-IR flux while 
representing little stellar mass. This includes, dusty stars 
of all luminosities and PAH emission. They do this by 
using independent component analysis (ICA) to sepa- 
rate out statistically significant source distributions and 
maximize [3.6] — [4.5] color differences. Using the 8.0 /jm 
images, they were able to obtain a rough estimate of 
the contribution of stellar and non-stellar (dust) sources. 
They note a selection bias against spatially coincident 
dust and clusters dominated by intermediate age stars 
and also a possible bias from mismatched [3.6] and [8.0] 
PSFs which would confuse the division between stellar 
and non-stellar emission. They caution the distinction 
between 3.6 /jm emission due to evolved stars and PAH 
emission are good for "rough estimates" of the compo- 
nents and that a more detailed analysis of the "contam- 
inates" is possible b ut beyond their scope. However, 
Meidt et al. (2012b) explores evolved stars in M100 in 
more detail. 

For comparison, [Melbourne et al.] ( |2012[ ) estimate the 
near-IR (HST WFC3/1R EI60W) H-bancf contributions 
of AGB stars in three of our galaxies, M81, NGC 300 and 
NGC 2403. The estimates are only for a single WFC3 
fields rather than global estimates. They find a NIR AGB 
luminosity contribution of 8% for M81, 15% for NGC 300 
and 17% for NGC 2403. Since the 1.6 /im emission is 
dominated by the photosphere of the AGB stars rather 
than the reprocessed emission from dust a round stars, 
these fractions need not match ours. Like Krick et al. 



(2010), they find that the stellar population models over- 
predict the AGB NIR luminosity and that the 1.6 /j,m 
AGB luminosity fraction increases with the mass fraction 
of intermediate age stars ( < 2 Gyr and < 0.3 Gyr) with 
no obvious metallicity trends. 

The "red stars" region is a mixture of sources. Most 
should be massive, evolved s tars such as 77 Carinae 
dHumphreys fc Davidson||1994[ ) IRC+10420 (e.g . |Jones 
et al.||l993|), or M33 Object X |Khan et al.||2011[ ) which 



are obscured by dust formed in mass ejections, borne are 
young star clusters containing dust, s uch as the br i ghtest 



3.6 /im source in M33 examined by Khan et al. (20111 



in the process of finding Object X. Exp losive transients 
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Fig. 5. — The fraction of the 3.6/im (filled triangles) and 
4.5/im (open circles) luminosity produced by stars in the red star 
region as a function of galaxy luminosity (top left), stellar mass 
(top right), specific star formation rate (lower left) and metallicity 
(lower right). 



Szczygiel et al. 2012). Since all these sources are likely 
massive stars or young clusters, we expected that the 
emissions from these red objects would be tied to the 
SFR. We see no correlation in our full sample between 
the bright, red source luminosity contribution and any 
galaxy property, including the SFR or SSFR. However, 
if we ignore the S MC, then there is a r elatively clean 



trend with SSFR. Bonanos et al. (2010) examined op 



can also form dust and lie in this region ( Fox et al.|2011 



tically selected luminous stars m the LMC and SMC, 
finding that those in the SMC seemed to be less dusty. 
Boyer et al. (2011 ) finds similar trends. Perhaps what we 
are observing is a combination of metallicity and SSFR, 
with the numbers rising with SFR but with dust produc- 
tion truncated at the lowest metallicity. This causes the 
bright, red stars of the SMC to be less red and dimmer 
in these IR bands. On average, the contribution of these 
sources is small, with 1.2% of the 3.6 /im flux and 1.9% 
of the 4.5 /im flux due to these bright, but rare obscured 
sources. 

5. DISCUSSION 

We surveyed the luminous dusty stars in 6 nearby 
galaxies. Using 3.6 and 4.5 /im data we identified AGB 
and other luminous red sources with significant emission 
from hot circumstellar dust. We determine their con- 
tribution to the flux at these bands and whether there 
is any correlation between their luminosity contribution 
and galaxy metallicity, stellar mass, luminosity or SSFR. 
Such studies can help to constrain one of the uncertain- 
ties considered in Conroy et al. (2010), namely the frac- 
tion of stellar emission obscured by stellar dust. By se- 
lecting the stars in the mid-IR we automatically include 
heavily obscured sources that would be missed in opti- 
cal studies. Additionally, by looking at whole galaxies, 
we have a relatively large sample of these relatively rare 
stars and should obtain reliable population statistics. 

The contributions of obscured AGB stars range from 
0.9% (M81) to 5.6% (SMC) at 3.6 /im and from 1.0% 
(M81) to 10.1% (SMC) at 4.5/im. We see a trend of 
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higher mid-IR with metallicity, SSFR, stellar mass and 
luminosity, but this is most likely a correlation with 
metallicity since metallicity is correlated with the other 
properties and none of the other properties should be 
correlated with AGB stars. The brighter, red sources 
contribute less, typically accounting for ~1.2% of the 
3.6 /Ltm and ~1.8% of the 4.5 /im luminosity. These 
sources are a mixture of young massive stars and star 
clusters, and we expect their luminosity fraction could 
be strongly correlated with the SSFR. The correlation 
may exist if there is also strong suppression of the red 



stars at the low met allicity of the SMC ( [Bonanos et al. 



2009 Bonanos et al.|[20T0 ). We lack a large enough dis- 



tribution of galaxy properties to address this question. 
These sources are probably more important at longer 
wavelengths. The 3.6 /im and 4.5 /im emissions are due 
to fairly hot dust close to the stars, while many of these 
sources are episodic dust sources rather than having long 
lived winds. During most of the period where they have 
a significant dust optical depth, the dust is at larger dis- 
tances and cooler temperatures radiating at longer wave- 
lengths. Simple estimates suggest that their contribu- 
tion will peak ^24 /mi, similar to 77 Carinae's present 



son 



day SEP (Robinson et al. 1973 Humphreys & David 



1994 1 . Dusty young star clusters will also te nd to 



be dominated by cooler dust ( Whelan et al.|[20lT |. Our 



findings are in agreement with recent results that while 
AGB stars have a significant effect on a galaxy's SED in 
the IR, the co ntribution is not as large as thought ini 



oT20 

Meidt et al ||2(Jl2b[ >. |Boyer et a'.| ( |20rT| 

stars contribute about 20% of the 3.6 



tially thought (|Kelson_fc Holden|20l"ol |Boyer et al.|2011 
' "" l find that AGB 



/im flux in the 
SMC, compared to our 6%. Our mid-IR identification 
of sources, likely selects a smaller, dustier sample of stars 
causing this difference. 

To put the impact of the IR light from these dusty 
stars into the broader context of the galaxy as a whole, 
we can compare our measurements to the bolometric lu- 
minositie s of the galaxies. Using the galaxy SED tem- 
plates of Assef et al. (2010) to estimate the bolometric 
luminosity of the galaxies, the summed 3.6 and 4.5 /im 
emission represents ~5% of the "bolometric" luminos- 
ity from (0.1-30 /im), so the summed 3.6 and 4.5 /mi 
emission of dusty AGB stars ranges from ~0.04% (M81) 
to ^0.37% (SMC) of the bolometric luminosity. If we 
model the mid-IR fluxes of the AGB stars as a stellar 
black body surrou nded by a dusty shell using DUSTY 
(Ivezic et al.| |1999p , we can estimate their total mid-IR 
luminosity. VVe find preferred dust temperatures of 1000- 
1500 K, as expected for emissions from dusty AGB winds 
with significant visual and negligible mid-infrared opti- 
cal depths. In these models, the summed 3.6 and 4.5 /tni 
AGB luminosity represents ~50% of the hot dust emis- 
sion from these stars. This means that the total dust 
emission of the AGB stars represents of order ~ 0.1% 
(M81) to - 0.8% (SMC) of the total luminosity of the 
galaxies. The difference between the observed 3.6 and 
4.5 /im emission and the total is modest because the dust 
is relatively hot. We cannot accurately estimate a similar 
correction for the red sources because they likely peak at 
a longer wavelengths (^24 /im) and we cannot effectively 
estimate a dust temperature. Overall, we find that cir- 
cumstellar dust around evolved stars represents a modest 
correction for galaxy SEDs but can be a significant cor- 



rection to their mid-IR colors. 

The uncertainty in the treatment of the AGB phase 
in the SPS models can have a si gnificant impact on the 



estimation of galaxy properties ( Marast on et al. 2006 
Conroy et al.|[2009{ Conroy et al.||2010| |(Jonroy Gunn 



2010 1. For example, Conroy et al.|(|2010[ ) find that lower 
metallicity populations prefer cooler AGB stars which 
could in part reflect increased dust production at lower 
metallicity. They also find that both TP- AGB stars and 
dust cause significant uncertainties in the properties o f 
star forming galaxies. In studies like |Conroy et al. ( 2010 ), 
it would be straight forward to include a reasonable phe- 
nomenological model of AGB or other stars with dusty 
winds. For a period T W i n d stars are assigned a dust op- 
tical depth r with the dust radiating at Tdust ~ 1000 K. 



The SEP c an easily be calculated using DUSTY (Ivezic 
| et al.|1999| or similar dust radiation transfer mod els. Ex- 
tending Ene UV t hrough near-IR wavelen gths of |Conroy 



et al. ( |2010[ ) and |Conroy fc Gunn| ( |20T0[ ) to the mid-IR 
would then constrain the hot dust contribution and its 
effects on the overall SED because it will be difficult for 
other variables to mimic the mid-IR dust emission. Fo- 
cusing on the 4.5 and 5.8 /im bands would minimize the 
need to worry about the strong PAH emissions at 8 /im 
and the weaker emissions at 3.6 /im. Although the overall 
amount of flux from dusty stars may be small, the effect 
on the mid-IR can be substantial and should be explored 
further. Our results are primarily limited by the low res- 
olution of Spitzer. The James Webb Space Telescope will 
allow this type of study to be performed on many more 
galaxies and over a broader wavelength range, provid- 
ing better constraints on how hot dust around stars, and 
AGB stars in particular, effect a galaxy's SED and how 
their contribution correlates with galaxy properties such 
as SFH and metallicity. 
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TABLE 1 
Galaxy Properties 



ID 


Distance 


Mctallicity 


Stellar Mass 


SFR 


SSRF 




Mpc 


12+log(0/H) 


Iog(M„/M ) 


Msyr- 1 


Gyr" 1 


M33 


0.96 


8.91 


8.79 


0.25 


0.41 


M81 


3.65 


8.99 


10.06 


0.46 


0.04 


NGC 300 


2.00 


8.65 


8.61 


0.12 


0.29 


NGC 2403 


3.22 


8.62 


9.00 


0.47 


0.47 


SMC 


0.06 


7.98 


7.29 


0.04 


1.54 


LMC 


0.05 


8.39 


8.00 


0.24 


1.96 



Note. — The distances are from Zaritsky ct al. (1994). expect 
for M33 jBonanos et al.|[2006t and M8TlfGerTre^aTjljll| l . The 
metallicities are i rom Zaritsky ct al. (1994) except tor th e SMC 
llPeimbert fc Torres- Permbert|1976j and LMC | |PageI et al.|1978| >. 
T he stellar masse s were determined using K band luminositi es Irom 
|Daleet al.| | |2009[ | and the mass to light ratio from |Bell et al.|(|2003| l. 
Vor the aMG and LM C the stellar masses were determined using 
the |DaIe et al.| (|2009[) [3.61 magnitude and an averaged mass to 
light ratio irom lalcon-Barroso et al.|(|2011j. The s tar formation 
rates wer e calculated using 11a me asurements from |Kennicutt et] 
[al] l |2008| with the conversion from |Meurer et al.| | |20TOf 



TABLE 2 
Total Galaxy Luminosity 



ID 


Raw (Lq) 


Corrected (Lq) 




[3.6] 


[4.5] 


[3.6] 


[4.5] 


M33 


3.78x10 s 


2.70x10 s 


(2.88 ±0.01) x 10 s 


(1.49 ± 0.004) x 10 s 


M81 


3.93x10 s 


1.93X10 9 


(3.84 ± 0.01) x 10 9 


(1.86 ±0.003) x 10 9 


NGC 300 


1.60 x10 s 


8.80X10 7 


(1.47 ±0.09) x 10 8 


(7.60 ±0.502) x 10 7 


NGC 2403 


5.59x10 s 


3.16x10 s 


(4.89 ±0.09) x 10 8 


(2.61 ± 0.041) x 10 s 


SMC 


2.81X10 7 


1.49X10 7 


(1.61 ± 0.17) x 10 7 


(8.09 ±1.124) x 10 6 


LMC 


1.11x10 s 


5.69X10 7 


(5.91 ± 0.23) x 10 7 


(2.98 ±0.131) x 10 7 



TABLE 3 
CMD Limits 



Region Limits 

MS slant: M 4 5 < 4.375([3.6] - [4.5]) - 10.563 

MS line: [3.6] - [4.5] < 0.10 

blue SM limit: [3.6] - [4.5] > -0.75 

MS upper limit: M4.5 > -13.60 

MS lower limit: M4.5 < -8.75 

red upper limits: M4 5 > —15.50 

red upper slant: M 4 . 5 > -2.727([3.6] - [4.5]) - 13.327 

AGB bottom slant: M4.5 < -1.271([3.6] - [4.5]) - 8.229 

AGB bottom limit: M 4 5 < -9.50 

AGB top slant: M4.5 > -3.437([3.6] - [4.5]) - 9.781 

AGB top limit: M 4 5 > -11.50 
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TABLE 4 

Background Contamination 



ID AGB Red Blue 





[3.6] 


[4.5] 


[3.6] 


[4.5] 


[3.6] 


[4.5] 


M33 


11 


11 


3 


2 


28 


31 


M81 


31 


31 


28 


26 


30 


32 


NGC 300 


20 


19 


11 


7 


48 


47 


NGC 2403 


43 


42 


29 


31 


31 


30 


SMC 


19 


18 


93 


84 


61 


61 


LMC 


32 


32 


21 


17 


62 


62 



Note. — These are the luminosities in the background regions 
as a percentage of the luminosity of the galaxy after correcting for 
area. 



TABLE 5 
Raw Luminosity Fractions 



ID 


AGB 


region 


Red stars 


Blue stars 




[3.6] 


[4.5] 


[3.6] 


[4.5] 


[3.6] 


[4.5] 


M33 


4.29 ±0.11 


7.37 ± 0.19 


1.42 ±0.16 


2.09 ±0.24 


14.96 ±0.39 


13.50 ±0.36 


M81 


0.86 ±0.01 


1.44 ± 0.02 


0.65 ± 0.03 


0.98 ±0.05 


2.12 ±0.04 


2.06 ± 0.04 


NGC 300 


3.59 ±0.12 


6.29 ±0.21 


1.06 ±0.17 


1.98 ±0.35 


7.78 ±0.35 


7.31 ±0.33 


NGC 2403 


4.27 ± 0.08 


6.88 ±0.13 


4.45 ±0.22 


6.20 ± 0.39 


8.59 ±0.27 


7.78 ± 0.25 


SMC 


7.49 ±0.57 


13.66 ± 1.10 


1.00 ± 0.47 


1.26 ±0.63 


55.35 ±3.39 


53.90 ±3.56 


LMC 


7.44 ±0.28 


13.86 ±0.53 


2.13 ±0.30 


3.24 ± 0.47 


44.09 ± 1.24 


41.24 ± 1.19 


Note. - 


The fractions are expressed 


as a percent of the total 






3.6 or 4.5/^m luminosities from Table [2] 
















TABLE 6 










Background-Corrected Luminosity Fractions 




ID 


AGB 


Red 


Blue 




[3.6] 


[4.5] 


[3.6] 


[4.5] 


[3.6] 


[4.5] 


M33 


2.85 ±0.10 


4.94 ± 0.17 


1.02 ±0.12 


1.54 ± 0.18 


8.00 ±0.72 


7.02 ±0.72 


M81 


0.58 ±0.02 


0.97 ±0.03 


0.46 ± 0.04 


0.71 ±0.06 


1.47 ± 0.09 


1.37 ± 0.09 


NGC 300 


2.71 ±0.15 


4.83 ±0.26 


0.89 ± 0.18 


1.73 ± 0.35 


3.80 ± 1.00 


3.67 ± 0.89 


NGC 2403 


2.38 ±0.14 


3.98 ±0.23 


3.07 ± 0.38 


4.27 ± 0.69 


5.81 ±0.54 


5.40 ± 0.46 


SMC 


5.55 ±0.86 


10.12 ± 1.81 


0.06 ±0.67 


0.18 ± 0.80 


19.61 ± 6.96 


19.06 ± 6.84 


LMC 


5.03 ±0.52 


9.47 ±0.98 


1.71 ± 0.44 


2.66 ±0.67 


16.63 ±3.35 


15.65 ± 3.07 



Note. — The fractions are expressed as a percent of the total 
3.6 or 4.5/^m luminosities from Table [2] 
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